Photonic Bloch oscillations and Zener tunneling in two-dimensional optical lattices 
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We report on the first experimental observation of pliotonic Blocli oscillations and Zener tunneling 
in two-dimensional periodic systems. We study the propagation of an optical beam in a square 
photonic lattice superimposed on a refractive index ramp, and demonstrate the tunneling of light 
from the first to the higher-order transmission bands of the lattice bandgap spectrum, associated 
with the spectral dynamics inside the first Brillouin zone and corresponding oscillations of the 
primary beam. 
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One of the fundamental concepts of quantum mechan- 
ics is a duality of particles and waves. Electromagnetic 
waves propagating in a periodic dielectric medium can 
behave in a rather similar fashion to electrons in a crys- 
talline potential. Hence many effects originally predicted 
in solid-state physics can be observed by monitoring light 
propagation in photonic structures. 

Two well-known fundamental phenomena associated 
with the propagation of waves and quantum particles in 
periodic media under the action of an external driving 
force are Bloch oscillations [2| and Zener tunneling U 
Indeed, in a periodic potential the dynamics of a par- 
ticle is dramatically affected by the bandgap structure 
of its energy spectrum. Moreover, due to the existence 
of the maximum wavenumber determined by the edge of 
the Brillouin zone, even particles associated with a sin- 
gle spectral band can behave in a rather unexpected way. 
For example, they do not follow the direction of a driv- 
ing force, but instead perform an oscillatory motion, the 
so-called Bloch oscillations. Such oscillations occur be- 
cause the presence of an external force causes the particle 
to gain momentum and approach the Bragg resonance. 
Since the band structure is periodic, the initial field dis- 
tribution is recovered after one crossing of the Brillouin 
zone. Therefore, one observes an oscillatory motion but 
no net shift of the particle in a real space. 

In order to explain why electrons move in a crystalline 
lattice under the action of a dc electric field, this single- 
band picture of the Bloch oscillations should be extended 
to take into account coupling to other bands, the effect 
known as Zener tunneling [3- If the variation of the 
superimposed linear potential within a unit cell is com- 
parable with the size of the gap to the adjacent band, an 
intra-band transition occurs. The strongest Zener tun- 
neling takes place when the particle reaches the edge of 
the first Brillouin zone where the gap is smallest. 

While originally predicted in the context of electrons 



in crystals, Bloch oscillations and Zener tunneling were 
also extensively investigated in different physical systems, 
including electrons in semiconductor superlattices 0, Q 
and cold atoms in optical lattices |^ . Recent progress 
in the fabrication and investigation of complex optical 
nanostructures has allowed for direct experimental ob- 
servations of many phenomena related to the wave prop- 
agation and interference. One-dimensional optical Bloch 
oscillations were observed in dielectric structures with a 
transversely superimposed linear ramp of the refractive 
index 0, @, S ■ A periodic distribution of the refrac- 
tive index plays a role of the crystalline potential, and 
the index gradient acts similar to an external force in 
a quantum system. It causes the beam to move across 
the structure where it experiences Bragg reflection on 
the high-index and total internal reflection on the low- 
index side of the structure, resulting in an optical ana- 
logue of Bloch oscillations. These experiments have re- 
cently been followed by the observations of tunneling in 
one-dimensional waveguide arrays [llLfl^ . It was demon- 
strated that for a small spectral gap even with a moderate 
strength of the index gradient, reasonably high transmis- 
sion rates from the first to the higher-order bands can be 
achieved facilitating a direct observation of the optical 
Zener tunneling. 

All direct experimental observations of Bloch oscilla- 
tions and Zener tunneling performed so far were limited 
to one-dimensional geometries. New effects may be asso- 
ciated with these phenomena in systems of higher dimen- 
sionality |l3l Il4j . In a two-dimensional periodic poten- 
tial the wave follows complex Lissajous-type trajectories 
when the direction of the static force does not coincide 
with a principal axis of the lattice. Additionally, the pro- 
cess of Zener tunneling becomes nontrivial as the band- 
gap structure can cause an enhanced tunneling in pre- 
ferred directions determined by the lattice symmetries. 

In this Letter, we report on the first experimental ob- 
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FIG. 1: (color online) (a) Interference pattern of four beams 
forming a square lattice; the inset shows symmetry points in 
the first Brillouin zone, (b) Experimentally measured profile 
of the background illumination Im (dash-dotted), and cor- 
responding change in the refractive index (solid). Dashed 
line indicates the initial position of the probe beam, (c) Cal- 
culated band structure of the lattice for Im = 0. (d) Band 
structure of the lattice vs. Im- Horizontal line is the adiabatic 
value of propagation constant; the dashed segment indicates 
Zener tunneling through the gap. 

servation of Bloch oscillations and Zener tunneling in 
two-dimensional (2D) periodic systems. By launching a 
laser beam into an optically-induced lattice with a super- 
imposed index gradient, we observe the occurrence of 2D 
Bloch oscillations and tunneling of light into the second 
transmission band of the lattice. 

The lattice is created by interfering four mutually co- 
herent ordinary-polarized broad beams in a biased pho- 
torefractive crystal 0, 0, 0| . The periodic light in- 
tensity distribution inside the crystal has the form of a 
square lattice shown in Fig. QJa) , described by the ex- 
pression 



Ig{x,y) = cos^ijrX/d) cos^iirY/d), 



(1) 



where {X, Y) = {x ^ y) / \/2, x and y are the transverse 
coordinates, A is the normalized amplitude, and d is the 
lattice period. Such a periodic light pattern induces a 
2D modulation of the refractive index for the extraordi- 
nary polarized probe beam. The ordinary polarized lat- 
tice beams, however, remain stationarity along the whole 
length of the crystal 0, 0| . Orientation of the lattice is 
set to 45° with respect to the c-axis of the crystal (hori- 
zontal in all figures; coincides with x axis in our notation) 
in order to reduce the effect of the intrinsic anisotropy of 
the photorefractive nonlinearity fl^. In order to create 
a transverse refractive index gradient we illuminate the 
crystal from the top with an incoherent white light which 
is modulated transversely, but is constant along the crys- 
tal length. The measured transverse intensity profile of 
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FIG. 2; (color online) (top) Numerically calculated position 
of the beam's center of mass; solid line - Bloch oscillations on 
the lattice; dashed line - beam deflection without the lattice. 
The vertical line indicates the position of the exit facet of the 
crystal. The light intensity distributions corresponding to the 
positions (a), (b), and (c) are depicted in the three snapshots. 

the white light is shown in Fig. ^b) (dash-dotted line) 
and it can be well approximated as 



/„(a;) =B[l + tanh(a;/77)]/2, 



(2) 



where the parameter r\ determines the extent of the in- 
duced index ramp. The total induced refractive index 
pattern is then given by [Fig.^b), solid] 



An(x,y) = 7 



Ig{x,y) +Im{x) 
1 + Ig(x,y) + Iyn{x) ' 



(3) 



where 7 can be tuned by varying the bias voltage. Then, 
we simulate the propagation of a probe optical beam by 
solving paraxial equation for the normalized electric field 
envelope E{x,y, z), 
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An[x,y)E^Q, (4) 



where z is the propagation distance. We use the following 
parameters to match our experimental conditions: wave- 
length in vacuum A = 0.532 /.tm, linear refractive index 
of the medium no = 2.35, lattice period d = 20/im and 
amplitude A — 0.5, 7 — 6.5, 77 — 100/im, and B = 1. 

Propagation of waves in a homogeneous lattice (for a 
fixed Im) is characterized by dispersion of Bloch waves 
which profiles have the form E = ^(X, Y) eyi\){iKxX + 
iKyY +il3z), where K = {Kx, Ky) is the Bloch wavevec- 
tor, the eigenvalue (3 is the propagation constant, and 
eigenmode profile ^' has the same periodicity as the un- 
derlying lattice. We show the first Brillouin zone in the 
inset of Fig. ^a), and the dispersion curves character- 
izing the band-gap spectrum are presented in Fig. ^c) 
for Im = 0, where the high-symmetry points are marked 
for the first two bands. The background illumination of 
the crystal, Im{x), results in the index distribution which 
represents a periodic structure imposed onto a monoton- 
ically increasing background. Because of the saturable 
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FIG. 3: (color online) Light intensity profiles (top row) and corresponding Fourier spectra (bottom) measured at the crystal 
output for different initial inclinations of the input beam. The inclination is (a) —0.095°, (b) —0.035°, (c) 0°, (d) 0.035°, and 
(d) 0.12°. The arrow in (a) indicates the direction of the index gradient. White diamonds indicate the first Brillouin zone. 



character of the photorefractive nonlinearity, the contrast 
of the resulting index grating varies transversely across 
the crystal, and the complete gap decreases, finally clos- 
ing for high background intensities, see Fig. ^d). Ex- 
periments were performed in the area of largest gradient 
where the gap is still open. 

We now model the propagation of an input Gaus- 
sian beam with a full width at half intensity maxi- 
mum of 37.5 /xm launched straight along the lattice, and 
plot in Fig. Eltop) the evolution of the beam center of 
mass, defined as / x\E\'^dr/ J \E\'^dr. As expected, the 
beam oscillates but, in contrast to the canonical Bloch- 
oscillations , the trajectories are deformed due to the 
transverse variations of the depth of the periodic poten- 
tial. According to the adiabatic theory, the Zener tunnel- 
ing from the initial Fi point in Figs.^c,d) occurs when 
the effective propagation constant (3 reaches the gap edge 
Ml and then tunnels through the gap to the point M2 , see 
the horizontal line in Fig. dd) . Due to the square sym- 
metry of the lattice this process results in the splitting of 
the initial beam into four parts, as shown in Figs.|2Ia-c). 
The intensity maxima of the strongest (reflected) por- 
tion are centered on the high-index points of the lattice. 
Hence this beam belongs to the fundamental band, and 
it keeps propagating in an oscillatory fashion. The other 
three beams are formed as a result of the tunneling to 
the second band, which is confirmed by the fact that 
their maxima are located in between the index maxima 
of the lattice. Weak or even vanishing gap on the high- 
index side promotes strong tunneling of light even for a 
moderate strength of the index gradient. 

We study experimentally the oscillations and tunnel- 
ing in 2D optical lattices using a 23 mm long, stron- 
tium barium niobate (SBN:65) photorefractive crystal, 
biased with an electric field of 5000 V/cm. A light from 
a green cw laser at 532 nm is used to create the lat- 
tice and the probe beam, while additional while light 
illuminator is used to induce the transverse index gradi- 



ent. First, we test the effect only of the induced index 
gradient on the beam propagation by switching off the 
lattice and launching the extraordinary polarized signal 
beam into the region with the steepest gradient. Mea- 
surements of the field at the output facet of the crystal 
performed with a CCD camera showed that the beam 
experiences a maximum transverse shift of 350 fim. In 
the next step, we launch the signal beam in the presence 
of optical lattice and the index gradient. Respective im- 
ages of the signal beam at the exit facet of the crystal are 
displayed in Figs. |3fa-e, top row). It is not possible to 
follow the evolution of the beam inside the crystal, but 
we can infer the details of its behavior by varying the 
incident angle of the probe beam. For angles below the 
Bragg resonance, the excitation by the probe with dif- 
ferent transverse wavenumber is equivalent to different 
starting points in the Brillouin zone. Therefore, when 
we launch the beam at different input angles we can scan 
the field evolution through the Bloch oscillation at the 
output facet of the crystal. In this way, we are able to 
monitor different stages of the Bloch oscillations, which 
would normally require crystals of different lengths. In 
Fig. I^Ja-e, top row) we can clearly see the predicted re- 
shaping of the beam and the tunneling of the beam en- 
ergy into three different channels. 

It should be noted that, as we tilt the input beam, we 
lose information about the absolute position of the Bloch 
trajectory. Nevertheless, we can study the inter-band 
tunneling by comparing the output intensity profiles. Ad- 
ditionally, fully conclusive picture of the beam evolution 
during a Bloch oscillation and the related Zener tunnel- 
ing can be obtained by monitoring the Fourier spectrum 
of the output field. This is realized by placing a lens in 
the beam path at the position of the image plane of the 
back side of the crystal. Then the focal plane of the lens 
contain the corresponding Fourier image. In this way 
various stages of the tunneling process can be resolved 
[see Fig. O bottom row]. The plot in Fig. |3Ia) corre- 



4 




FIG. 4: (color online) Light intensity profiles calculated (a-c) 
and measured (d-f ) at the crystal output for a single site exci- 
tation. (a,d) discrete diffraction with no transverse gradient, 
(b,e) weak and (c,f) strong index gradients. 

sponds to the initial beam tilt of —0.095° (negative an- 
gle corresponds to an initial motion of the beam against 
the force produced by the gradient). For this angle the 
beam intensity profile at the output facet of the crystal is 
strongly modulated with adjacent maxima being approx- 
imately out-of-phase. Corresponding intensity distribu- 
tion in the Fourier domain shows that the beam propa- 
gates at the Bragg angle (0.32°) in the M-symmetry point 
of the lattice. Increasing the input angle allows us to scan 
through different positions inside the Bloch period and to 
monitor the light tunneling into the higher-order bands. 
Evidently, only one peak lies inside the first Brillouin 
zone for all the incident angles, and it corresponds to 
the field which undergoes Bloch oscillations. The other 
three peaks representing tunneled beams lie outside the 
boundaries of the first Brillouin zone, thus they belong to 
the second transmission band. For an angle of —0.035° 
[Fig. Ol^b)] the beam just experiences its first Bragg re- 
flection and the three tunneled beams emerge from all 
M-symmetry points of the lattice. At larger angles all 
the beams get accelerated again by the index gradient, 
and the central part of the beam completes a full Bloch 
oscillation, see Figs. Ol^c-e). We note that changing the 
angle between the lattice and the index gradient leads 
to asymmetry in scattered and tunneled beams, but the 
qualitative picture remains the same. 

As discussed above, the Zener tunneling accompany- 
ing Bloch oscillations of a broad beam that covers several 
sites of the lattice occurs only when the beam approaches 
the edge of the Brillouin zone. A drastically different sce- 
nario has been predicted for an input beam which is 
initially very narrow, comparable in its size with a single 
site of the 2D lattice. Such a beam excites simultane- 
ously waves with wavevectors distributed over the whole 
Brillouin zone. This results in a symmetric breathing of 
the beam in the first band as it periodically diffracts and 
refocuses in propagation 0. Because the narrow beam 
excites modes with the wave vectors distributed over the 



whole Brillouin zone, there always exists a component 
which propagates in the vicinity of the band edge. Hence, 
the tunneling occurs now continuously over the whole 
Bloch period. To observe this effect we use a tightly fo- 
cused probe beam (20/im) which is centered at one of 
the lattice sites and launched into the crystal at normal 
incidence. We record the light intensity distribution at 
the exit facet of the crystal varying the strength of the 
transverse index gradient. Figures ^a,d) show the dis- 
crete diffraction of a narrow beam in the lattice with no 
superimposed index gradient. Note the almost perfect 
spatial symmetry of the field distribution which reflects 
the geometry of the lattice. In the case of a rather weak 
gradient [Figs.^b,e)] the light travels approximately two 
thirds of the Bloch period before reaching the end of the 
crystal. The central part of the beam starts to refocus 
while some light already tunnels to the second band. Al- 
though no pronounced spatial separation of the light in 
the first and second bands occurs at the available crystal 
length, the induced asymmetry of the field distribution 
is a clear indication that tunneling takes place. For a 
strong gradient in Figs.^c,f) almost no breathing is vis- 
ible and the light in the first band is confined to the input 
waveguide and its direct neighbors. 

In conclusion, we have studied the propagation of light 
beams in a two-dimensional optically-induced photonic 
lattice in the presence of a transverse ramp of the refrac- 
tive index. We have observed two-dimensional photonic 
Bloch oscillations and Zener tunneling to distinct points 
of the bandgap spectrum. 

This work was produced with the assistance of the Aus- 
tralian Research Council under the ARC Centres of Ex- 
cellence program. 
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